ValueFolderTemp='/sps/isc/rneveu/EtudeIRM/Nature/ResultsJournal/';
VecValueTime=[0 2 1 5 10 3 4];
VecValueFilter=[0 1];
VecResultsJournal=[];
VecResultsJournal2=[];
for n1=1:33
    load([ValueFolderTemp 'FileAnalysisCharacteristicsEditorsResearchJournal' num2str(n1) '.mat']);
    VecResultsJournal=[VecResultsJournal;VecResultsJournalEditor];
end
for n5=1:size(VecValueTime,2)
    for n6=1:size(VecValueFilter,2)
        ValueTime=VecValueTime(1,n5);
        ValueFilter=VecValueFilter(1,n6);
        n3=0;
        n6=0;
        VecResults=[];
        VecResults4=[];
        VecResults3=[];
        VecResults5=[];
        VecResultsCorrelation5=[];
        VecResults52=[];
        VecResults7=[];
        VecResultsJournalTemp=[];
        VecResultsJournal2Temp=[];
        VecCorrelationJournal=NaN(21,41);
        VecCorrelationJournal2=NaN(21,41);
        VecValueFigure=[];
        VecResults6={};
        VecResults6(1,:)={'Parameter' 'mean (SD)' 'n (%)' 'Cohen effect size d' 'p-value ttest' 'p-value Wilcoxon test' 'mean (SD) before' 'mean (SD) after'};
        VecResults6(4,1)={'Slope value (robustfit) of N original articles with a date in the journal/N total articles with a date in the journal'};
        VecResults6(5,1)={'Slope value (regress) of N original articles with a date in the journal/N total articles with a date in the journal'};
        VecResults6(2,1)={'Intercept (robustfit) of N original articles with a date in the journal/N total articles with a date in the journal'};
        VecResults6(3,1)={'Intercept (regress) of N original articles with a date in the journal/N total articles with a date in the journal'};
        VecResults6(10,1)={'Slope value (robustfit) of N articles with a date in the journal with no authors with a previous publication in all Nature journals/N original articles with a date in the journal'};
        VecResults6(11,1)={'Slope value (regress) of N articles with a date in the journal with no authors with a previous publication in all Nature journals/N original articles with a date in the journal'};
        VecResults6(12,1)={'Slope value (robustfit) of N articles with a date in the journal with no authors with a previous publication in the journal/N original articles with a date in the journal'};
        VecResults6(13,1)={'Slope value (regress) of N articles with a date in the journal with no authors with a previous publication in the journal/N original articles with a date in the journal'};
        VecResults6(6,1)={'Intercept (robustfit) of N articles with a date in the journal with no authors with a previous publication in all Nature journals/N original articles with a date in the journal'};
        VecResults6(7,1)={'Intercept (regress) of N articles with a date in the journal with no authors with a previous publication in all Nature journals/N original articles with a date in the journal'};
        VecResults6(8,1)={'Intercept (robustfit) of N articles with a date in the journal with no authors with a previous publication in the journal/N original articles with a date in the journal'};
        VecResults6(9,1)={'Intercept (regress) of N articles with a date in the journal with no authors with a previous publication in the journal/N original articles with a date in the journal'};
        VecResults6(18,1)={'Slope value (robustfit) of N articles with a date in the journal with at least one intermediate author with a previous publication in all Nature journals/N original articles with a date in the journal'};
        VecResults6(19,1)={'Slope value (regress) of N articles with a date in the journal with at least one intermediate author with a previous publication in all Nature journals/N original articles with a date in the journal'};
        VecResults6(20,1)={'Slope value (robustfit) of N articles with a date in the journal with at least one intermediate author with a previous publication in the journal/N original articles with a date in the journal'};
        VecResults6(21,1)={'Slope value (regress) of N articles with a date in the journal with at least one intermediate author with a previous publication in the journal/N original articles with a date in the journal'};
        VecResults6(14,1)={'Intercept (robustfit) of N articles with a date in the journal with at least one intermediate author with a previous publication in all Nature journals/N original articles with a date in the journal'};
        VecResults6(15,1)={'Intercept (regress) of N articles with a date in the journal with at least one intermediate author with a previous publication in all Nature journals/N original articles with a date in the journal'};
        VecResults6(16,1)={'Intercept (robustfit) of N articles with a date in the journal with at least one intermediate author with a previous publication in the journal/N original articles with a date in the journal'};
        VecResults6(17,1)={'Intercept (regress) of N articles with a date in the journal with at least one intermediate author with a previous publication in the journal/N original articles with a date in the journal'};
        VecResults6(22,1)={'Difference in number of original articles published in Nature and Nature Communications the year before the launch of a new journal and the year after the launch'};
        VecResults6(23,1)={'Difference in number of original articles published in Nature the year before the launch of a new journal and the year after the launch'};
        VecResults6(24,1)={'Difference in number of original articles published in Nature and Nature Communications over the two years before the launch of a new journal and over the two years after the launch'};
        VecResults6(25,1)={'Difference in number of original articles published in Nature over the two years before the launch of a new journal and over the two years after the launch'};
        VecResults6(26,1)={'Number of original articles published in the new journal the year of its launch'};
        VecResults6(27,1)={'Number of original articles published in the new journal over the two years after its launch'};
        VecResultsCorrelation(1,:)={'Parameter' 'r pearson, mean (SD)' 'p-value ttest' 'p-value Wilcoxon test' '% of variance explained' 'r spearman, mean(SD)' 'p-value ttest' 'p-value Wilcoxon test' 'n (%)'};
        VecResultsCorrelation(2,1)={'Correlation between rate of original articles with authors with no previous article published in all Nature journals and rate of original articles'};
        VecResultsCorrelation(3,1)={'Correlation between rate of original articles with authors with no previous article published in the journal and rate of original articles'};
        VecResultsCorrelation(4,1)={'Correlation between rate of original articles with at least one intermediate author with an article published previously in all Nature journals and rate of original articles'};
        VecResultsCorrelation(5,1)={'Correlation between rate of original articles with at least one intermediate author with an article published previously in the journal and rate of original articles'};
        VecResultsCorrelation(6,1)={'Correlation between number of original articles and total number of articles'};
%         n1=32;
%         if ValueTime>=1
%             if ValueFilter==0
%                 ValueFile=['FileAnalysis' num2str(ValueTime) 'Journal' num2str(n1) '.mat'];
%             else
%                 ValueFile=['FileAnalysis' num2str(ValueTime) 'Filtered' num2str(ValueFilter) 'Journal' num2str(n1) '.mat'];
%             end
%         else
%             if ValueFilter==0
%                 ValueFile=['FileAnalysisJournal' num2str(n1) '.mat'];
%             else
%                 ValueFile=['FileAnalysisFiltered' num2str(ValueFilter) 'Journal' num2str(n1) '.mat'];
%             end
%         end
%         load([ValueFolderTemp ValueFile]);
%         n2=1;
%         TestTemp=0;
%         while n2<=size(VecResultsArticle2,1)&&TestTemp==0
%             if VecResultsArticle2(n2,4)>0
%                 TestTemp=1;
%             end
%             n2=n2+1;
%         end
%         VecValueTemp32=VecResultsArticle2(n2-1:end,:);
%         n1=33;
%         if ValueTime>=1
%             if ValueFilter==0
%                 ValueFile=['FileAnalysis' num2str(ValueTime) 'Journal' num2str(n1) '.mat'];
%             else
%                 ValueFile=['FileAnalysis' num2str(ValueTime) 'Filtered' num2str(ValueFilter) 'Journal' num2str(n1) '.mat'];
%             end
%         else
%             if ValueFilter==0
%                 ValueFile=['FileAnalysisJournal' num2str(n1) '.mat'];
%             else
%                 ValueFile=['FileAnalysisFiltered' num2str(ValueFilter) 'Journal' num2str(n1) '.mat'];
%             end
%         end
%         load([ValueFolderTemp ValueFile]);
%         n2=1;
%         TestTemp=0;
%         while n2<=size(VecResultsArticle2,1)&&TestTemp==0
%             if VecResultsArticle2(n2,4)>0
%                 TestTemp=1;
%             end
%             n2=n2+1;
%         end
%         VecValueTemp33=VecResultsArticle2(n2-1:end,:);
        for n1=1:33
            if ValueTime>=1
                if ValueFilter==0
                    ValueFile=['FileAnalysis' num2str(ValueTime) 'Journal' num2str(n1) '.mat'];
                else
                    ValueFile=['FileAnalysis' num2str(ValueTime) 'Filtered' num2str(ValueFilter) 'Journal' num2str(n1) '.mat'];
                end
            else
                if ValueFilter==0
                    ValueFile=['FileAnalysisJournal' num2str(n1) '.mat'];
                else
                    ValueFile=['FileAnalysisFiltered' num2str(ValueFilter) 'Journal' num2str(n1) '.mat'];
                end
            end
            ValueFileTemp=dir([ValueFolderTemp ValueFile]);
            if size(ValueFileTemp,1)==1
                load([ValueFolderTemp ValueFile]);
                n2=1;
                TestTemp=0;
                while n2<=size(VecResultsArticle2,1)&&TestTemp==0
                    if VecResultsArticle2(n2,4)>0
                        TestTemp=1;
                    end
                    n2=n2+1;
                end
                VecValueTemp=VecResultsArticle2(n2-1:end,:);
%                 if n1<32
%                     if size(VecValueTemp,1)>size(VecValueTemp33,1)-1
%                         if size(VecValueTemp,1)>1
%                             VecResults7=[VecResults7;VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1),4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1,4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1),4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1,4) sum(VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)-1:size(VecValueTemp32,1)-size(VecValueTemp,1),4)) sum(VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1:size(VecValueTemp32,1)-size(VecValueTemp,1)+2,4)) sum(VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)-1:size(VecValueTemp32,1)-size(VecValueTemp,1),4)) sum(VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1:size(VecValueTemp32,1)-size(VecValueTemp,1)+2,4)) VecValueTemp(1,1) sum(VecValueTemp(1:2,1)) n1];
%                         else
%                             VecResults7=[VecResults7;VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1),4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1,4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1),4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1,4) NaN NaN NaN NaN VecValueTemp(1,1) NaN n1];
%                         end
%                     else
%                         if size(VecValueTemp,1)>1
%                             VecResults7=[VecResults7;VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1),4)+VecValueTemp33(size(VecValueTemp33,1)-size(VecValueTemp,1),4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1,4)+VecValueTemp33(size(VecValueTemp33,1)-size(VecValueTemp,1)+1,4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1),4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1,4) sum(VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)-1:size(VecValueTemp32,1)-size(VecValueTemp,1),4)+VecValueTemp33(size(VecValueTemp33,1)-size(VecValueTemp,1)-1:size(VecValueTemp33,1)-size(VecValueTemp,1),4)) sum(VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1:size(VecValueTemp32,1)-size(VecValueTemp,1)+2,4)+VecValueTemp33(size(VecValueTemp33,1)-size(VecValueTemp,1)+1:size(VecValueTemp33,1)-size(VecValueTemp,1)+2,4)) sum(VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)-1:size(VecValueTemp32,1)-size(VecValueTemp,1),4)) sum(VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1:size(VecValueTemp32,1)-size(VecValueTemp,1)+2,4)) VecValueTemp(1,1) sum(VecValueTemp(1:2,1)) n1];
%                         else
%                             VecResults7=[VecResults7;VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1),4)+VecValueTemp33(size(VecValueTemp33,1)-size(VecValueTemp,1),4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1,4)+VecValueTemp33(size(VecValueTemp33,1)-size(VecValueTemp,1)+1,4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1),4) VecValueTemp32(size(VecValueTemp32,1)-size(VecValueTemp,1)+1,4) NaN NaN NaN NaN VecValueTemp(1,1) NaN n1];
%                         end
%                     end
%                 end
                VecValue2Temp=VecResultsArticle1(n2-1:end,:);
                VecResults2=VecValueTemp;
                VecValueTemp2=1:size(VecValueTemp,1);
                VecValueTemp2=VecValueTemp2';
                VecValueTemp2=VecValueTemp2-1;
                if size(find(VecValueTemp(:,1)>0),1)>0
                    VecResults4=[VecResults4 VecResultsArticle2(:,4)];
                end
                if size(VecValueTemp2,1)>2
                    VecValueTemp5=[VecValueTemp(:,3)./VecValueTemp(:,4)*100 VecValueTemp(:,1)./VecValueTemp(:,3)*100 VecValueTemp(:,2)./VecValueTemp(:,3)*100];
                    VecValueTemp51=[VecValue2Temp(:,3)./VecValue2Temp(:,4)*100 VecValue2Temp(:,1)./VecValue2Temp(:,3)*100 VecValue2Temp(:,2)./VecValue2Temp(:,3)*100];
                    ValueTemp5=find(VecValueTemp(:,4)==0);
                    VecValueTemp5(ValueTemp5,:)=zeros(size(ValueTemp5,1),3);
                    ValueTemp5=find(VecValueTemp(:,3)==0);
                    VecValueTemp5(ValueTemp5,:)=zeros(size(ValueTemp5,1),3);
                    ValueTemp5=find(VecValue2Temp(:,4)==0);
                    VecValueTemp51(ValueTemp5,:)=zeros(size(ValueTemp5,1),3);
                    ValueTemp5=find(VecValue2Temp(:,3)==0);
                    VecValueTemp51(ValueTemp5,:)=zeros(size(ValueTemp5,1),3);
                    ValueResults=robustfit(VecValueTemp2,VecValueTemp5(:,1));
                    ValueResults2=regress(VecValueTemp5(:,1),[ones(size(VecValueTemp2,1),1) VecValueTemp2]);
                    ValueResults3=robustfit(VecValueTemp2,VecValueTemp5(:,2));
                    ValueResults4=regress(VecValueTemp5(:,2),[ones(size(VecValueTemp2,1),1) VecValueTemp2]);
                    ValueResults5=robustfit(VecValueTemp2,VecValueTemp5(:,3));
                    ValueResults6=regress(VecValueTemp5(:,3),[ones(size(VecValueTemp2,1),1) VecValueTemp2]);
                    ValueResults31=robustfit(VecValueTemp2,VecValueTemp51(:,2));
                    ValueResults41=regress(VecValueTemp51(:,2),[ones(size(VecValueTemp2,1),1) VecValueTemp2]);
                    ValueResults51=robustfit(VecValueTemp2,VecValueTemp51(:,3));
                    ValueResults61=regress(VecValueTemp51(:,3),[ones(size(VecValueTemp2,1),1) VecValueTemp2]);
                    ValueCorrelation1=corr(VecValueTemp5(:,1),VecValueTemp5(:,2),'type','pearson');
                    ValueCorrelation2=corr(VecValueTemp5(:,1),VecValueTemp5(:,2),'type','spearman');
                    ValueCorrelation3=corr(VecValueTemp51(:,1),VecValueTemp51(:,2),'type','pearson');
                    ValueCorrelation4=corr(VecValueTemp51(:,1),VecValueTemp51(:,2),'type','spearman');
                    ValueCorrelation5=corr(VecValueTemp5(:,1),VecValueTemp5(:,3),'type','pearson');
                    ValueCorrelation6=corr(VecValueTemp5(:,1),VecValueTemp5(:,3),'type','spearman');
                    ValueCorrelation7=corr(VecValueTemp51(:,1),VecValueTemp51(:,3),'type','pearson');
                    ValueCorrelation8=corr(VecValueTemp51(:,1),VecValueTemp51(:,3),'type','spearman');
        %             ValueResults31=robustfit(VecValueTemp(:,4),VecValueTemp(:,3));
        %             ValueResults31=regress(VecValueTemp(:,3),[ones(size(VecValueTemp,1),1) VecValueTemp(:,4)]);
                    ValueCorrelation9=corr(VecValueTemp(:,4),VecValueTemp(:,3),'type','pearson');
                    ValueCorrelation10=corr(VecValueTemp(:,4),VecValueTemp(:,3),'type','spearman');
%                     if isnan(ValueResults5)==1
%                         if sum(VecValueTemp(:,3))==0
%                             ValueResults5=0;
%                             ValueResults6=0;
%                         end
%                     end
                    VecResults5=[VecResults5;ValueResults(1,1) ValueResults2(1,1) ValueResults(2,1) ValueResults2(2,1) ValueResults3(1,1) ValueResults4(1,1) ValueResults31(1,1) ValueResults41(1,1) ValueResults3(2,1) ValueResults4(2,1) ValueResults31(2,1) ValueResults41(2,1) ValueResults5(1,1) ValueResults6(1,1) ValueResults51(1,1) ValueResults61(1,1) ValueResults5(2,1) ValueResults6(2,1) ValueResults51(2,1) ValueResults61(2,1) n1];
                    VecValueTimeTemp=1:size(VecValueTemp5,1);
                    VecValueTimeTemp=VecValueTimeTemp';
                    VecValueFigure=[VecValueFigure;VecValueTemp5 n1*ones(size(VecValueTemp5,1),1) VecValueTemp51 VecValueTimeTemp repmat(VecResultsJournal(n1,:),size(VecValueTemp5,1),1)];
                    VecResultsJournal2Temp=[VecResultsJournal2Temp;VecValueTemp5(end,:) VecResultsJournal(n1,3) VecResultsJournal(n1,16) VecResultsJournal(n1,27)];
                    VecResultsCorrelation5=[VecResultsCorrelation5;ValueCorrelation1 ValueCorrelation2 ValueCorrelation3 ValueCorrelation4 ValueCorrelation5 ValueCorrelation6 ValueCorrelation7 ValueCorrelation8 ValueCorrelation9 ValueCorrelation10];
                    VecResultsJournalTemp=[VecResultsJournalTemp;n1];
%                     figure
%                     plot(VecValueTemp2,VecValueTemp5(:,1),'-')
%                     hold on
%                     plot(VecValueTemp2,VecValueTemp5(:,2),'-r')
            %         hold on
            %         plot(VecValueTemp2,VecValueTemp5(:,3),'-c')
            %         figure
            %         plot(VecValueTemp2,VecValueTemp(:,3)/max(VecValueTemp(:,3)),'-')
            %         hold on
            %         plot(VecValueTemp2,VecValueTemp(:,1)/max(VecValueTemp(:,1)),'-r')
            %         figure
            %         plot(VecValueTemp2,VecValueTemp(:,3),'-')
            %         hold on
            %         plot(VecValueTemp2,VecValueTemp(:,4),'-r')
                end
                n3=n3+1;
                n6=n6+1;
            else
                VecResults52=[VecResults52;n1];
            end
        end
        VecValueTemp4=[zeros(1,2) ones(1,2) zeros(1,4) ones(1,4) zeros(1,4) ones(1,4)];
        for n1=2:21
            VecResults6(n1,2)={[num2str(mean(VecResults5(:,n1-1)),3) ' (' num2str(sqrt(var(VecResults5(:,n1-1))),3) ')']};
            VecResults6(n1,3)={[num2str(size(VecResults5,1),3) ' (' num2str(size(VecResults5,1)/n6*100,3) ')']};
            if VecValueTemp4(1,n1-1)==1
                VecResults6(n1,4)={num2str(mean(VecResults5(:,n1-1))/sqrt(var(VecResults5(:,n1-1))),3)};
                [ValueTemp,ValueTemp2]=ttest(VecResults5(:,n1-1));
                ValueTemp3=signrank(VecResults5(:,n1-1));
                VecResults6(n1,5)={num2str(ValueTemp2,3)};
                VecResults6(n1,6)={num2str(ValueTemp3,3)};
            else
                VecResults6(n1,4)={'-'};
                VecResults6(n1,5)={'-'};
                VecResults6(n1,6)={'-'};
            end
            VecResults6(n1,7)={'-'};
            VecResults6(n1,8)={'-'};
        end
%         for n1=1:4
%             ValueTemp2=find(isnan(VecResults7(:,(n1-1)*2+1))==0);
%             VecResults6(n1+21,2)={[num2str(mean(VecResults7(ValueTemp2,(n1-1)*2+2)-VecResults7(ValueTemp2,(n1-1)*2+1)),3) ' (' num2str(sqrt(var(VecResults7(ValueTemp2,(n1-1)*2+2)-VecResults7(ValueTemp2,(n1-1)*2+1))),3) ')']};
%             VecResults6(n1+21,3)={[num2str(size(ValueTemp2,1),3) ' (' num2str(size(ValueTemp2,1)/31*100,3) ')']};
%             VecResults6(n1+21,4)={num2str(mean(VecResults7(ValueTemp2,(n1-1)*2+2)-VecResults7(ValueTemp2,(n1-1)*2+1))/sqrt(var(VecResults7(ValueTemp2,(n1-1)*2+2)-VecResults7(ValueTemp2,(n1-1)*2+1))),3)};
%             [ValueTemp,ValueTemp2]=ttest(VecResults7(ValueTemp2,(n1-1)*2+2)-VecResults7(ValueTemp2,(n1-1)*2+1));
%             ValueTemp3=signrank(VecResults7(ValueTemp2,(n1-1)*2+2)-VecResults7(ValueTemp2,(n1-1)*2+1));
%             VecResults6(n1+21,5)={num2str(ValueTemp2,3)};
%             VecResults6(n1+21,6)={num2str(ValueTemp3,3)};
%             VecResults6(n1+21,7)={[num2str(mean(VecResults7(ValueTemp2,(n1-1)*2+1)),3) ' (' num2str(sqrt(var(VecResults7(ValueTemp2,(n1-1)*2+1))),3) ')']};
%             VecResults6(n1+21,8)={[num2str(mean(VecResults7(ValueTemp2,(n1-1)*2+2)),3) ' (' num2str(sqrt(var(VecResults7(ValueTemp2,(n1-1)*2+2))),3) ')']};
%         end
        for n1=1:5
            VecResultsCorrelation(n1+1,2)={[num2str(mean(VecResultsCorrelation5(:,(n1-1)*2+1)),3) ' (' num2str(sqrt(var(VecResultsCorrelation5(:,(n1-1)*2+1))),3) ')']};
            [ValueTemp,ValueTemp2]=ttest(VecResultsCorrelation5(:,(n1-1)*2+1));
            ValueTemp3=signrank(VecResultsCorrelation5(:,(n1-1)*2+1));
            VecResultsCorrelation(n1+1,3)={num2str(ValueTemp2,3)};
            VecResultsCorrelation(n1+1,4)={num2str(ValueTemp3,3)};
            VecResultsCorrelation(n1+1,5)={[num2str(mean(VecResultsCorrelation5(:,(n1-1)*2+1).^2),3) ' (' num2str(sqrt(var(VecResultsCorrelation5(:,(n1-1)*2+1).^2)),3) ')']};
            VecResultsCorrelation(n1+1,6)={[num2str(mean(VecResultsCorrelation5(:,(n1-1)*2+2)),3) ' (' num2str(sqrt(var(VecResultsCorrelation5(:,(n1-1)*2+2))),3) ')']};
            [ValueTemp,ValueTemp2]=ttest(VecResultsCorrelation5(:,(n1-1)*2+2));
            ValueTemp3=signrank(VecResultsCorrelation5(:,(n1-1)*2+2));
            VecResultsCorrelation(n1+1,7)={num2str(ValueTemp2,3)};
            VecResultsCorrelation(n1+1,8)={num2str(ValueTemp3,3)};
            VecResultsCorrelation(n1+1,9)={[num2str(size(VecResultsCorrelation5,1),3) ' (' num2str(size(VecResultsCorrelation5,1)/33*100,3) ')']};
        end
%         [ValueTemp,stats]=robustfit(VecResultsJournal2Temp(:,4),VecResultsJournal2Temp(:,2));
%         [ValueTemp2,stats2]=robustfit(VecResultsJournal2Temp(:,5),VecResultsJournal2Temp(:,2));
%         [ValueTemp3,stats3]=robustfit(VecResultsJournal2Temp(:,6),VecResultsJournal2Temp(:,2));
        [ValueTemp,ValueTemp2]=corr(VecResultsJournal2Temp(:,2),VecResultsJournal2Temp(:,4),'type','spearman');
        [ValueTemp3,ValueTemp4]=corr(VecResultsJournal2Temp(:,2),VecResultsJournal2Temp(:,5),'type','spearman');
        [ValueTemp5,ValueTemp6]=corr(VecResultsJournal2Temp(:,2),VecResultsJournal2Temp(:,6),'type','spearman');
%         VecResultsJournal2=[VecResultsJournal2;ValueTemp' stats.p' NaN ValueTemp2' stats2.p' NaN ValueTemp3' stats3.p' NaN ValueTime ValueFilter];
        VecResultsJournal2=[VecResultsJournal2;ValueTemp ValueTemp2 ValueTemp3 ValueTemp4 ValueTemp5 ValueTemp6 ValueTime ValueFilter];
%         save(['/sps/isc/rneveu/EtudeIRM/Nature/AnalysisThirdLevelJournalTime' num2str(ValueTime) 'Filtered' num2str(ValueFilter) '.mat'],'VecResults6','VecResultsCorrelation');
        save([ValueFolderTemp 'AnalyseMixedModelJournalTime' num2str(ValueTime) 'Filtered' num2str(ValueFilter) 'Research.mat'],'VecValueFigure');
    end
end
        